Spectroscopy of capacitively coupled Josephson-junction qubits 
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We show that two capacitively-coupled Josephson junctions, in the quantum limit, form a simple 
coupled qubit system with effective coupling controlled by the junction bias currents. We compute 
numerically the energy levels and wave functions for the system, and show how these may be tuned 
to make optimal qubits. The dependence of the energy levels on the parameters can be measured 
spectroscopically, providing an important experimental test for the presence of entangled multiqubit 
states in Josephson-junction based circuits. 
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Ramos et al. have proposed that electrically well- 
isolated Josephson junctions can be used as qubits [ij. 
Two recent experiments using different isolation schemes 
have reported Rabi oscillations in single junctions 0, 
demonstrating the existence of macroscopic quantum 
coherence. While longer coherence times are desir- 
able, these experiments show that single Josephson junc- 
tions are strong candidates for solid-state qubits; several 
Josephson-based types have been proposed Q . 

One of the next major steps towards building a 
Josephson-junction based quantum computer will be the 
observation of quantum properties of coupled macro- 
scopic qubits. A simple scheme for making coupled 
qubits, junctions connected by capacitors, has recentTy 
been proposed by Blais et al. and Ramos et al. 
This scheme is illustrated for the two-qubit case in 
Fig. 1(a). 

In this paper, we focus on the immediately accessible 
fundamental experiments-spectroscopic measurements of 
macroscopically entangled quantum states-that are pos- 
sible with this system. We calculate, using highly accu- 
rate numerical methods, the energy levels and metastable 
wave functions for the circuit of Fig. 1(a) in terms of the 
junction parameters, bias currents, and coupling capac- 
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FIG. 1: (a) (left) Circuit diagram for two idealized capac- 
itively coupled Josephson junctions, (b) (right) The tilted 
washboard potential for a single current-biased Josephson 
junction with three metastable quantum states. 



itance. Our numerical analysis demonstrates that the 
system can be tuned to create appropriately spaced en- 
ergy levels and coupled states. The features that we dis- 
cuss can both guide the experimental effort of observing 
multiqubit quantum states and provide help in optimiz- 
ing the design of qubits and gates. But we emphatically 
stress that experimental observation of these macroscopic 
entangled quantum states will be an important achieve- 
ment in its own right, and will provide strong support 
for the validity of macroscopic quantum mechanics and 
the existence of macroscopic entanglement l6|. Spectro- 
scopic observation of these states should be possible by 
using standard single-junction experimental techniques 
iff. 

The Hamiltonian for an ideal single current-biased 
Josephson junction, with critical current Ic and junction 
capacitance Cj, is 

7J(7,p) - (AEch-^) p^-Ej {cos J + Jj), (1) 

where 7 is the gauge-invariant phase difference across 
the junction, J = I /Ic, I is the (tunable) bias current, 
Ej — (<I>o/c/27r) is the Josephson energy, Ec — e^/2Cj is 
the charging energy, and $0 — h/2ei& the flux quantum. 
The canonical momentum isp = (<I>o/27r) Cj7 = hQ/2e^ 
where Q is the charge on the junction. The ratio Ej /Ec 
determines whether the system is in a phase, charge, or 
intermediate regime. The qubits explored in this paper 
have Ej » Ec and hence are in the phase regime. 

The equations of motion for a single current-biased 
junction are equivalent to those for a particle in the 
tilted washboard potential shown in Fig. 1(b). Classi- 
cally, for J < 1 there are stable minima about which the 
phase can oscillate with the characteristic plasma fre- 
quency COp (J) = v/27r/c/$oCj (1 " J^)^^" 0. Quantum 
mechanically, the system exhibits localized metastable 
states in each well that can tunnel out into the run- 
ning (finite-voltage) state. The effective barrier height 
^Ubarrier [scc Fig. 1(b)] for a siuglc junction in units of 
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FIG. 2: (a) Potential V' with strong coupling (C = 0.8) and 
Ji = J2. The coupling induces a squeezing in the 7+ direction, 
relative to the 7- direction, (b) V' contours with ( = 0.8 and 
Ji 7^ J2, around the vicinity of one well. The symmetry of 
V' shows that, despite detuned bias currents, 7+ and 7_ are 
approximately the normal modes, (c) V' with small coupling 
((" — 0.01) and Ji 7^ J2 showing (by symmetry) that 71 and 
72 are approximate normal modes, (d) V' with ^ — 0.01 and 
Ji = J2 showing that 7+ and 7_ are normal modes. 



htOp (J) , assuming J < 1, is related to Ns by 
A[/, '^3/4 / 1/2 



^barrier 



hujp (J) 



(2) 



Here, iVg is the approximate number of metastable bound 
states for a single isolated junction Our analysis ex- 
plores the relevant regime for quantum computing where 
Ns is small and the nonlinearity of the Hamiltonian is 
important. 

By adjusting the bias current it is possible to tune the 
barrier height to obtain, for example, three metastable 
energy levels Eq < Ei < E2 with the two lowest 
states forming the basis |0) , |1) of a qubit. State |2) 
has the highest escape rate due to tunneling and can 
therefore act as an auxiliary readout state, where read- 
out is achieved by microwave pumping at a frequency 
UJ12 = {E2 — El) /h. Detection of a voltage across the 
junction implies that the system was previously in the 
state |1) and has entered the running state. 

The Hamiltonian for the coupled two-junction circuit 
shown in Fig. 1(a) is 



H 



AEc 



{pI+pI + 2CpiP2) 



(1 + 0^2 
Ej (cos 71 -I- cos 72 -I- Ji7i -I- J272) , 



(3) 



where C = Cc'/{Cc + Cj) is the dimensionless cou- 
pling parameter, Cc is the coupling capacitance, and 
J1.2 are the normalized bias currents of junctions 1 
and 2, respectively. The canonical momenta, pi 2 — 
(Cc + Cj) ($o/27r)^ (71,2 - C72,i) , are proportional to 
the charges on each junction plus the charge on the cou- 
pling capacitor plate adjacent to it 0,0- 

By making a canonical change of variables, defined by 



7± 



(7i±72)/V2(l±C), 



p± = V2(l±C)bi±P2) 



(4) 
(5) 



find the transformed Hamiltonian 



H' (p+,_p-,7+,7-) = 



^Ec 



(p^+/_)+^'(7+,7-). 

(6) 

Here, the momentum coupling term 2(^pip2 in the orig- 
inal Hamiltonian is shifted to coupling in the new po- 
tential energy V — y(74.,7_). Figures 2(a,b) illustrate 
how the coupling induces a squeezing of V along the 7+ 
direction; a strong coupling of ^ = 0.8 has been chosen 
to accentuate this behavior. 

We gain further insight into the coupling of the junc- 
tion states by looking at the normal modes for small os- 
cillations. For small coupling and detuned bias currents 
( Ji far from J2) the normal modes are approximately 71 
and 72, i.e., the junctions are effectively decoupled. This 
effect is shown in Fig. 2(c) by the approximate symme- 
try of V with respect to reflections about the 71 and 72 
axes. When Ji = J 2 the normal modes become 7+ and 
7_, and we therefore expect the coupled junction states 
to be entangled symmetric and antisymmetric combina- 
tions of the single junction states. This can be seen in 
Fig. 2(d), where V is symmetric with respect to reflec- 
tions about the 7+ and 7_ axes. Figure 2(b) shows V 
for C, = 0.8 and unequal bias currents; despite detuning 
the large C, prevents the junctions from decoupling and 
7+ and 7-- effectively remain as normal modes. 

The challenging demands of quantum computing re- 
quire an accurate and precise quantitative description of 
the states going beyond simple perturbation theory. To 
achieve this, we have computed the states and energy 
levels numerically using a nonperturbative fast Fourier 
transform split-operator method ^ applied to the full 
nonlinear Hamiltonian in Eq. (jj)). Our implementation 
computes the wave functions on a lattice using a fourth- 
order integration of the imaginary-time evolution opera- 
tor exp(— ffr). While this method is relatively slow, its 
results for a subset of system parameters conflrm that the 
much faster complex scaling method applied to the 
cubic approximation of the full potential Iy| is accurate 
to at least 0.1%. The faster complex scaling method then 
allows us to compute energy levels for a wide range of sys- 
tem parameters. A further important property of both 
these numerical methods is that they are well suited to 
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FIG. 3: The modulus squared of the probability amplitude 
of the first six quasistationary wave functions for capacitively 
coupled current-biased Josephson junctions with coupling of 
^ = 0.01 and normalized bias currents of 3\ = 3i= 0.986 93. 



finding metastable states in potentials that allow tunnel- 
ing, particularly in more than one dimension, where other 
methods fail. The computed quantum states have been 
further verified by time evolving them on a lattice using 
real time split-operator methods. This has confirmed dy- 
namically that these states are truly quasistationary and 
thus accurately determined; where applicable, agreement 
has also been found with higher-order WKB analysis. 

In Fig. 3 we show, for example, the numerically com- 
puted wave functions for identical junctions with capac- 
itances Cj — 4.3 pF and critical currents Ic — 13.3 /iA. 
Junctions with these physical characteristics are readily 
fabricated and of physical interest. Figure 3 shows the 
modulus squared of the wavefunctions of the first six qua- 
sistationary states for the coupling strength C — 0.01, 
and with bias currents J\ = Ji = 0.986 93 such that 
isolated junctions would have approximately three qua- 
sistationary states (Ng — 3). These large bias currents 
make the nonlinearities of the potential pronounced, and 
the states deviate significantly from coupled harmonic- 
oscillator states. The states \n) in Fig. 3 are ordered by 
energy a rounded bracket has been used to distin- 
guish the coupled two-junction states \n) from single- 
junction states \n). The second and third states ex- 
pressed in terms of single-junction direct product states 
are |1) (|01) - |10)) /V2 and |2) = (|01) + |10)) |^/2, 
whereas the higher-energy states are more complicated 
superpositions that depend upon the bias currents and 
coupling. The ordering of the states in Fig. 3 may 
be understood by looking at the potentials shown in 
Fig. 2(a,b); wave functions extended in the 7+ direction 
have higher energy because of the coupling induced ef- 
fective squeezing in the 7+ direction, relative to the 7_ 
direction. Observe that because the (71,72) configura- 
tion space variables are the collective degrees of freedom 
of distinct junctions, the wave functions represent macro- 
scopic nonlocally correlated (and hence entangled) states. 

For designing qubits out of coupled junctions we need 
to know how the energy levels depend on coupling and 
bias current. Figure 4 shows the effects of varying the 



111) 

I20> 






102) 


\T\ ■ 




110) 
101) 


|2)S (|10>+|01»/V2 




|1) = (I01>-|10»/V2 




ICQ) 


= 10) 









18 



0.05 0.1 0.15 0.2 

C = Cc/(Cc+C,) 



FIG. 4: Frequencies versus coupling strength for equal bias 
currents Ji = J2 = 0.986 93, Ic = 13.3 /xA, and Cj = 4.3 pF. 
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FIG. 5: Frequencies of the first six states versus bias current 
Ji with J2 = 0.986 93 fixed and a coupling strength oi C, = 
0.01, Ic = 13.3 and Cj = 4.3 pF. 

coupling strength in the range < C < 0.2 on the 
first six energy levels with Ji — J2 = 0.986 93. The 
plasma frequency of each single junction when C, — Q 
is LOp {Ji = J2) /27r = 6.2037 GHz. The states are la- 
beled at the left of Fig. 4 for zero coupling, where the 
product representation \nm) = \n) ® \m) is appropriate. 
For zero coupling the nonlinearity of the potential has 
broken the degeneracy between |5) = |11) and the pair 
(|3) = |02),|4) = |20)). 

In real experiments, the coupling strength C, will typ- 
ically be fixed by the circuit design, and can only be 
varied by making a completely new sample. By contrast, 
the bias currents through each junction are easily varied, 
and allow manipulation of the entangled states shown 
in Fig. 3. Figure 5 shows how the energy levels change 
for C = 0.01 and J2 = 0.986 93 fixed, while Ji is varied 
around J2. There are prominent avoided level crossings 
indicated in the figure for both on-tune (Ji = J2) and 
off-tune (Ji ^ J2) bias currents. The predicted gap for 
the |1) , |2) on-tune splitting at Ji = 0.986 93 is 57 MHz 
for Ic = 13.3 and Cj — 4.3 pF. The predicted gap 
for the first off-tune splitting between states |4) and |5) 
at Ji = 0.98444 is 80 MHz, and for the second off-tune 
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FIG. 6: Frequencies versus bias current Ji with J2 = 0.98 
fixed and a coupling strength of = 0.05, Ic = 13.3 fiA, and 
Cj = 4.3 pF. 

splitting at Ji = 0.989 62 the gap is 72 MHz. The pre- 
dicted gap for the on-tune |3) , |4) splitting of 4 MHz is 
much smaller than the others because it is a second-order 
avoided crossing in perturbation theory. Figure 6 shows 
the energy levels for the same junction parameters but 
with C = 0.05 and J2 = 0.98 (N^ ~ 5). 

We have labeled the states in Figs. 5 and 6 at the far 
left and right -when the currents are detuned and hence 
the states are effectively uncoupled~as product states. 
This labeling is only strictly correct when ^ = 0. The 
mixing that occurs between states when the bias currents 
are brought into tune is indicated for the states |1) and 
|2) in Figs. 4 and 5. A swaplike gate operation can be 
constructed by exploiting this mixing p|- 

Experimental data similar to Figs. 5 and 6 would be 
important first evidence for the existence of macroscopic 
entangled states like those shown in Fig. 3. A typical 
experiment to probe the energy levels in Figs. 4-6 would 
proceed by preparing the system in the ground state |0) 
by cooling well below T ~ E^i/k ~ 300 mK, where 
Enni = {En - Em) , and E„ = En{Ji, C) are the coupled- 
junction energies whose dependence on bias current and 
coupling we have shown in Figs. 4-6. Varying the bias 
current Ji (with J2 fixed) while simultaneously inject- 
ing microwaves at a frequency u) should lead to an en- 
hancement in the tunneling from the zero-voltage state 
to the finite-voltage running state of the system when 
^ = Enm (Ji) /fi- This enhancement produces a corre- 
sponding peak in escape rate measurements Q • By 
varying u) and Ji for the coupled junctions, we may map 
out the energy levels for comparison with Figs. 5 and 6. 

Experimentally, the expected energy gap between the 
avoided levels can be resolved if both the quality factor 
Q and coupling ^ of the system are reasonably large. For 
example, with the typical junction parameters assumed 
here, and with ( = 0.01, the predicted gap for the |1),|2) 
splitting of 57 MHz can easily be resolved with a Q of 200 
pl]. Furthermore, one can reasonably track the bending 



of the resonant escape peaks near the avoided crossings 
within the span of a typical experimental current window 
of about 30 nA 0. These splittings can be made even 
easier to detect by increasing the coupling capacitance 
(see Fig. 6), though large coupling may inhibit efficient 
quantum gates. 

In conclusion, we have presented predictions for fun- 
damental experiments that probe macroscopic entangled 
states by the relatively simple scheme of doing spec- 
troscopy on coupled junctions while varying external bias 
currents. The energy levels of these entangled states 
should be readily observable using the same experimental 
techniques that have allowed spectroscopy of single junc- 
tions. The numerical methods we have used are powerful 
tools for mapping out the metastable states of nonlinear, 
many-level coupled systems, and allow us to explore a 
wide range of junction parameters and couplings. This 
kind of detailed study will be necessary for the design of 
realistic coupled qubits. 
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